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We review the application of Migdal's Theory of Finite Fermi System to 
the structure of deformed nuclei, approaches beyond the conventional linear 
response, and microscopic calculations of the Migdal-parameters. 

1 Personal recollection by Josef Speth 

In 1975 I met Arkadi Benediktowitsch for the first time in Gerry Brown's institute in 
Stony Brook. He was interested in our numerical results of his theory of finite Fermi 
systems. During a visit in Dubna in 1977 he invited me to his apartment in Moscow 
for lunch and discussions with some of his collaborators. He visited our institute in 
Jiilich 1980 for a month, where he impressed not only physicists but also non- physicists 
at various parties in Jiilich. The photograph shows him as a magician who makes some 
conjuring trick for my children. 

2 Introduction 

There exist several reviews on Migdal's Theory Finite Fermi System jT] (TFFS) where 
the application of the theory to various aspects of atomic nuclei has been presented. The 
original approach is restricted to the solution of the linear response function, which is 
connected with the 4-point many-body Green function (GF). In the early review by Speth 
et al. P] numerical results for nuclei in the lead region have been summarized. They 
also reported on an extension to the 6-point GF, that had been previously developed 
in Ref. [21 !3- Khodel and Saperstein investigated in their publications, which have 
been summarized in Ref. [5] , a self-consistent version of the TFFS and presented results 
for several closed shell nuclei. Drozdz et al. [6] developed an important extension 
of the TFFS which allows to describe giant resonances in a quantitative way. The 
authors include Iplh and 2p2h configurations in a consistent manner and demonstrate 




Figure 1: A.B. Migdal as magician in Speth's iiouse (1980). 

the power of their method for electric, magnetic and charge exchange resonances. In a 
more recent work by Kamerdzhiev et al. [7] another extension of the TFFS was reviewed. 
Here Ipl/iCgiphonon configurations were included in the conventional linear response 
equation, their effect on collective electric resonances discussed and results for closed 
shell nuclei presented. Finally Tselyaev [51 E] presented a model where he included not 
only two phonon states but also pairing correlations in a consistent manner, which allows 
applications to open shell nuclei. With the exception of Ref. [H El E] mostly electric 
properties of spherical nuclei have been discussed. Recent applications to neutron rich 
nuclei are reviewed in Ref. [TOlIIIlIiaiiaElIISlIiaill]. 

In the present review we demonstrate that the TFFS has been successfully applied 
to strongly deformed nuclei where electric as well as magnetic properties in the rare 
earth and actinide region have been calculated. It is very important to mention that the 
derivation of the basic formulas within the many-body GF formalism is very general and 
the corresponding equations are not restricted to the TFFS. We will see in the next sec- 
tion that for the numerical application one needs single particle energies, single particle 
wave functions and an effective (ph)-interaction. If pairing is included one also needs an 
effective particle-particle (pp) interaction. In the TFFS the single particle quantities are 
deduced from phenomenological single-particle models and the residual interactions are 
parametrized. A comparison between the phenomenologically determined parameters 
and calculated ones from first principles can be found in Ref. [IB] . In the self consistent 
approaches the single particle data are obtained from a mean field calculated with an 
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effective Hamiltonian, Lagrangian or energy functional. The residual ph-interaction is 
given in such approaches as second derivative of the energy. In this review we want to 
show that the TFFS is not only able to reproduce known data but that most of the 
predictions have been experimentally verified. The review may also give a guide line for 
those scientists which start to investigate deformed nuclei, though in different models 
and theories, respectively. 

In section 2 we define the GF and the corresponding equations of motion. These 
quantities are derived with the functional method developed by Kadanoff and Baym [TU] 
and its generalization to systems with pairing proposed by Brenig and Wagner [20]. In 
the derivation of the basic formulas we do not give the details but we outline the basic 
idea and refer to other reviews or the original publications for details. In this section 
we also define the second-order response function which has been successfully used to 
calculate isomer shifts of rotational states in the rare earth region [21]. In section 3 
we present results for deformed nuclei. This involves electric as well as magnetic low- 
lying and high- lying resonances. The microscopic results are compared with the well 
known phenomenological models. We discuss in sections 4 an extension of Migdal's 
linear response theory, the second order response. This formalism was application to 
a very specific nuclear structure effect, namely the isomer shifts of rotational states in 
deformed nuclei. These extremely small effects were measured with the help of the 
Mossbauer effect. In section 4 we review magnetic states in deformed nuclei. In section 
5 we discuss microscopic calculations of the Migdal parameters from first principles and 
in section 6 a short summery is given. We refer mostly to the work done in Munich and 
Jiilich, with additional references to the work by Urin |22|I23]. who investigated collective 
vibrations in deformed nuclei within the TFFS in a different way and the application of 
the TFFS to the Inglis cranking model by Birbrair [24] . 

3 Method 

3.1 Many Body Green Functions 

The GF from which we start are defined as the ground-state expectation value of a 
N-particle system of the time order product of pairs of quasi-particle creation and an- 
nihilation operators. In the present context, these Green functions are functions of an 
external field g(l,2), which is a formal device to derive higher GF and their equations 
of motion [TO]. After one has performed the derivations one can put the source field 
equal to zero and arrives at the conventional definition. In a single particle basis the 
one-particle GF is given as: 



5,(1,2) 
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JJ — gij rf5d6 g(5,6)a+ (i5)a^g(i6) 



(2) 
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and the short hand notation: 

Y.ldti = I dl. (3) 

The two-particle and three-particle GF have correspondingly the form (with q(l,2)=0): 

5(13,24) ={^){A{}\Ta,,{ti)a,,{h)a+^{t2)aX{U)\A^) (4) 

5(135,246) = (i3) (A0| ra,,(ti)a,3(t3)aK5(*5)ai(t2)a+ (t4)a+ (te) 1^0) . (5) 

The equation of motion for the one-particle GF is the Dyson equation which has the 
form: 



j d3{5(l,3) + g(l,3) + S,(l,3)}g,(3,2) =5(1,2). (6) 



with the abbreviations: 

5u,.Ati-t2)[i^^-el}^ =5(1,2). (7) 

The quantity S is the self-energy or mass operator, an effective one-body potential, 
which in principle is given by the bare interaction of the corresponding Hamiltonian 
and the two- particle GF [2]. We define the linear response function L by the functional 
derivative due to the field (?(1, 2) of the one-particle GF: 

= L(13,24) = 5(13,24) - 5(1, 2)^(3, 4). (8) 

The functional derivative of the Dyson equation gives an integral equation for the re- 
sponse function 

L(13, 24) = -5(1,4)5(3,2) - i J (iSdedTdS 5(1, 5)K(57,68)L(83, 74)5(6, 2), (9) 
where we introduced the effective two-body interaction K via 

||i|^«(13,24), (10) 

The change of the one-particle GF 5gq due to an external field 6q is given in linear 
response : 

6gg{l,2) = J d3d4L{13,24)6q{4,3). (11) 

Analogously to the linear response function one defines the second order response func- 
tion as: 

A f Vfi = ^135,246) = 5(135, 246) + 25(1, 2)5(3, 4)5(5, 6) (12) 
-5(13, 24)5(5, 6) - 5(15, 26)5(3, 4) - 5(35, 46)5(1, 3). 
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As we will see in the following Eq.Q is the central equation in the conventional TFFS. 
The one-particle GF and the two-particle GF are defined self-consistently by a system of 
non-linear equations. This is, however, of little use for practical applications. In order to 
arrive at solvable equations, one applies Landau's quasi particle concept and his renor- 
malization procedure, which he developed for his Theory of Fermi Liquid [2F| . Migdal 
introduced this concept into the nuclear many-body problem. Here, the quasi particles 
are the single-particle states. Following Landau, one splits the (Fourier transformed) 
one-particle GF into a pole part and a remainder. Written in the configuration space of 
the single-particle wave functions '^y the equation has the form: 



9u^U2{() = , ^''^".^ I V + al^uii'^)- (13) 

eu^-e + irj s\gn[ey^ - /i) 

The are the single-particle energies, Zy the single particle strength and /j. is the Fermi 
energy. The main goal is to obtain an equation for the response function that can be 



solve in practice. With the ansatz in Eq. (13) one writes the product of two GF as a 
singular part S and a rest B: 

g{e,n)g{e) = S{e,n) + B{e,n). (14) 

The singular part has the form: 



here the rii^ are the occupation numbers for quasi particles: 1 and for particles and 
holes, respectively. 

3.2 Linear response 



From the linear response equation ( 9|11 ) on can derive an equation for the change of the 
density 5p due to an external field 6q: 

{ey, - - ^) ^Pv^v2 = {nv^ - ny^) (Sq^.y^i^) + ^ F^tiv^vz hv-iVi)- (16) 
Here PP^ and 5q are the renormalized p/i-interaction and renormalized external field. 



respectively, which includes the Zu factors as well as the regular part B of Eq. ( 14 ) . The 



homogeneous part of Eq.(16) is formally identical with the conventional Random Phase 
Approximation (RPA) 



From Eq.(17) one calculates excitation energy $7^ and the corresponding (renormalized) 
ph-transition amplitude Xj^i^a from the ground-state to the excited state \Am) of an 
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even-even nucleus with mass number A. From the latter ones follow the expectation 
values of one-body operators: 

{Am\Q\A0) = Y^Ql(l^xT2u,. (18) 

The renormalized single-particle operators Q'^^^ also include the z^, factors and the reg- 
ular part B. In the case of electric multipole operators Q*^-^-^ correspond to the bare 
operators due to charge conservation, whereas the magnetic operators are parametrized, 
with universal parameters 151. The response function includes also an equation for 
moments and transitions in the neighboring odd-mass nuclei [X^, 

{A±l,a\Q\A±l,f3) - <^,,;3 (^0| Q |ylO) = T^p{€ap,Q), (19) 

the corresponding equation for the vertex operator 
Ta/3{Gai3,Q) has the form: 

TuiU2{^al3,Q) = Qa/ SuiaSu2l3 + ^ Fu^UiU2U3 ^ T^-iuA^^af) , Q) ■ (20) 

u-xUA ^1^3 



Here ta^ is the energy difference between the two states a and /3. From Eq.(19) we 
see that in the case of moments one can only calculate the difference between the even 
and odd mass nuclei. This, however, allows a very precise calculation of the differences 
of charge distributions (isotope shifts) . 



3.3 Many Body Green functions including pairing correlations 

With the exception of closed shell nuclei, pairing correlations play an important role in 
nuclei. Therefore one needs an extension of the previously discussed Theory of Finite 
Fermi systems which includes pairing correlations. Such an extension has first been 
presented by Larkin and Migdal . Here we give a more general derivation [27\ |251 122] 
which is based on generalized GFs [30] . We introduce a one-particle GF matrix of the 
form: 



Qkl 



/ ^11 ^12 



r<2i r'22 



{A\T{U}\A) 



(21) 



{A\T{Ua^at}\A) - {A\T{Ua,aj} \A + 2) 

{A + 2| T{Ua±at} \A) -{A + 2\ T{Ua±aj} \A + 2) 

The bar on the indices denotes the time- reversed states. The two-particle GFs and the 
response functions are defined in the same way as before and are obtained as functional 
derivatives. The generalized Dyson equation possesses also four components which in- 
cludes four mass operators. By functional derivation of the Dyson equation one obtains 
the integral equations for the four response functions. The ansatz for the (diagonal) pole 
parts of the four one-particle GFs can be written in a compact form [21] : 
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gW(c.) = - ( ^ + p ] (22) 



L,= r" -"^^n ■,T,= ) (23) 

\^ uxvx -ui J ' ^ \ -uxvx -vx J 

with the BCS quantities: 

■"^-K'-^) "-50^^) ^^-(--">'-^'v (^^> 

The gap A is given by the usual gap equation: 

here F^^ is the renormahzed particle-particle (pp) interaction, which also enters in the 
equation for the response function. 

3.4 Quasi-particle RPA 

The four coupled equations for the response functions can be reduced to two coupled 



equations. With the ansatz for the pole part of the one-particle GFs given in Eq.(22) 
plus a regular part one performs an analog renormalization procedure as described in 
the previous section. The final equations have the form of the well known quasi-particle 
RPA (QRPA) equations which allow to calculate e.g. collective excitations in super fluid 
Fermi systems of two quasi-particle type. These equations have been previously derived 
in different ways by Bogoliubov Baranger [32) and Belyaev |3H]. Birbrair fI7\ and 
Kamerdzhiev |28j also used the GF formalism making explicit the difference between 
pp- and ph- interaction. Here we give the compact form of the equations derived by 
Baranger [32]: 

(^A + E^) Z+ + Y: {vLFf^U^^ HtF^;U^.) V = nZ^. (26) 

{Ex + E,) Zl^ + ^ (r/,-,<-, V Z+ = ^Zl^ (27) 

with the normalization condition: 

2E^A.(^A.)* = 1- (28) 
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The various quantities in Eqs.(26) and (27) are defined as follows: 



U\Un =F VxVk,] 



UXVk =F vxu,^, 



(29) 



and 



:pph± 



(30) 



From the solutions of Eqs. (26) and (27) we obtain the excitation energies O and 



the amplitudes which are connected with the transition probabilities B{EL) in the 
following way. 

2 



B {EL) = e2 (2 - 5ko) 



with 



K - Vxk^Xk^ Xxk - 2 V-^^K ^Ak ) ' 



(31) 



(32) 



where is the renormalized version of the following (unrenormalized) matrix element 
of an A-particle system: 



= (^0|a+a,|ylm) 



(33) 



Eq. (31) is the analog to Eq. (18) for super fluid systems. Like in the non super fluid 



case, the electric multipole operators can be replaced by the bare operators, because of 
charge conservation. In the magnetic case one has to use renormalized operators. The 
analog to Eq. (20) (moments and transitions in odd mass nuclei) is given in Ref. |34j- 
This equation has been successfully applied to the calculation of isomer shift [22] and 
isotope shifts [33] in deformed odd mass nuclei. 



3.5 Second order response theory 

Within the linear response theory one calculates the change of the expectation value 
of a single particle operator in an external field. Due to the linear relation, the single 
particle operators and the external field have to have the same operator structure. One 
of the especially nice applications of the following extended version of the TFFS is the 
calculation of the change of the nuclear charge radii in rotational states. These tiny 
effects were measured with help of the Mossbauer effects. Here one calculates within the 
cranking model the change of the scalar operators 5rp^ due to the Coriolis perturbation 
5q = — Jx which is a pseudo vector operator. Here ilc is the so called cranking 
parameter. 

The change of the density in linear response p^^^ due to the Coriolis perturbation has 
the form |21] : 
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)\2 



J-)A-) 



^ A/t. 



(1) 



(34) 



It is obvious that the change of the radius is zero in hnear response. This equation has 
been previously derived by Migdal [35] and Birbrair [36] and has been used to calculate 
moments of inertia and gyromagnetic ratios PT. The equation for the change of the 
density in second order p^^^ has the same structure as Eq. (34). All quantities that 



enter in the equation are know from the linear response theory with the exception of an 
effective three particle interaction which has been neglected in all applications. 



(2) ~(2) r. -,1 



(35) 



The inhomogeneous term depends in a complicated way quadratically on the Coriolis 
perturbation. The very lengthy formula for /)'^2^[inh] and the modified F'^^ and F'^^ are 
given in Ref. |21j . 



4 Application to deformed nuclei 

In order to solve the basic equations one needs as input single particle- wave functions, 
single-particle energies and the ph-and pp-interaction. Migdal has designed his theory 
in close connection to Landau's Fermi liquid theory. Therefore one takes as far as 
possible the input data from experiment or from models which reproduce the needed 
experimental data as good as possible. If the RPA equation is derived within the many- 
body Green function formalism [2j one obtains an explicit form for the ph-interaction 
F^^ that depends on the single-particle model and the size of the configuration space. 
In addition it is nonlocal and energy-dependent. In practice the ph-interaction is not 
calculated from that formula but it is parametrized. Following Landau's procedure F^^ 
is transformed into momentum space and one considers the interaction on the Fermi 
surface of nuclear matter. Here one can replace the energies by the Fermi energy and 
the magnitude of the momenta by the Fermi momentum. In this approximation F^^ 
depends only on the angle between the ph-momenta P and P' before and after the 
collision; (we suppress the spin and isospin dependence) 




Here Pi{x) is the Legendre polynomial of order / and the constants Fi are the famous 
Landau-Migdal parameters. One introduces dimensionless parameters by defining: 

Fi = Cofi. (37) 
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Here Cq is the inverse density of states at the Fermi surface. One choses density depen- 
dent parameters which correct for the finite size of the nuclei. For deformed nuclei an 
axially deformed Fermi distribution was used. 

• 0"2Tl • T2] (38) 

The pp-interaction is treated in the same way. In all applications so far, the interaction 
is restricted to the scalar-isoscalar component, with parameters which are in some cases 
density dependent. 

The single-particle wave functions are taken from a single-particle model and the 
single-particle energies are taken as far as possible from experiment. The deformed 
rare earth nuclei and the actinides are well described within the unified model [38 that 
provides a good working single-particle model. Here the authors obtained the single 
particle wave functions from a deformed Woods-Saxon potential [35]. The most critical 
question has been the choice of the single particle energies. In the case of the isomer 
shifts, which was the first application of the TFFS to deformed nuclei, the level scheme by 
Ref. ^D] has been used with some corrections due to new experimental informations. For 
the solution of the QRPA equation which had been done some years later the theoretical 
level schemes have been thoroughly readjusted to new experimental data. |41] . 

5 Results 

5.1 Classical quadrupole shape oscillations 

The intrinsic wave functions in deformed nuclei with axial symetry are characterized 
by the parity and the K-quantum number (projection of the angular momentum on 
the rotational axis). The low-lying K'^= 0"*" and K'^= 2+ has been investigated in the 
past in great detail. In the framework of the incompressible liquid drop model with 
axially symmetric equilibrium, these low-lying collective states have been interpreted as 
quadrupole shape oscillations. The K'^= 0^ vibrations which preseve axial symmetry are 
called /3-vibrations and the K'^= 2+ modes which break the axial symmetry are referred 
to as 7- vibrations. In addition one obtains also K'^= 1+ states which correspond to the 
collective rotation about an axis perpendicular to the symmetry axis. These states are 
usually called the spurious ones, as they are not connected with internal excitations. In 
microscopic calculations one obtains as many l"*" solutions as one has ph components. 
Only the lowest solution which is the most collective one is "spurious". All the other ones 
are internal excitations and correspond to real physical states as we will see in section 4.3 
section. The major part of the non-spurious (isoscalar) strength is concentrated around 
11 MeV. A schematic representation of the quadrupole vibration of the phenomenological 
liquid drop model are shown in Fig. ([2|. 

In section 4.4 the classical picture is compared with the microscopically calculated 
transition densities. 
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COLLECTIVE STATES 




K = i 2,Y-vibration 



Figure 2: Schematic representation of the classical quadrupole shape oscillations in a 
spheroidal nucleus. 

5.2 Low-lying electric states in the rare earth region 

Experimentally, collective K'"= 0^ excitations with energies of about 1 MeV, the so called 
/?- vibrations, have been known for a long time. The corresponding K'"= 2'^ excitations, 
the 7-vibrations, are also present in this energy range. The considerable variations of the 
energies and transition probabilities over the region of deformed nuclei may cast some 
doubt on the classical interpretation. Actually, the corresponding high-lying collective 
states, as we will see, correspond much more that interpretation. The low-lying states 
are dominated by only a few ph-components and depend therefore sensitively on the 
single-particle level scheme and the transition densities have little similarity with the 
classical picture. The theoretical energies and transition probabilities, which can be 
found in Ref. [41j, are in general in fair agreement with the data. 

5.3 High-lying electric states in the rare earth region 

The high-lying collective states are the well known giant resonances which are qualita- 
tively different in deformed nuclei compared to spherical nuclei. The phenomenological 
model predicts for a given multipolarity a splitting of the different K-components. This 
is experimentally well established for the electric dipole resonances t42j, where as for 
the quadrupole resonances a broadening is observed [13] . From QRPA calculations one 
obtains in a natural way the excitation energies, transition probabilities and the mag- 
nitude of the splitting between different K'^ components. As an example the results 
for the 2+ and 0+ states in ^™Yb are shown in Fig. ([3|. The calculation has been 
performed in a large, but discrete basis, therefore the theoretical states do not have any 
width. The B{EX) strength is summed in intervals of 0.5 MeV. An interesting results 
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Figure 3: Distribution of B{E2) and B{EO) strength in ^^Yb. The B{EX) strength is 
summed in intervals of 0.5 MeV. The peaks are identified by capital letters. A: 
lowest K'^ = 0+ excitation. B: lowest = 2+ excitation. C,D, and E: isoscalar 
giant quadrupole resonances for K'^ = 0+, 1+ and 2+ components, respectively. 
F: = 0+ state (AT = 0), predominately of a breathing mode type (the 
states C and F are both superpositions of /3- vibrations and breathing mode). 
C,D and F: isovector quadrupole resonances. J: This K"^ = 0+ state (AT=1 ) is 
predominantly an isovector breathing mode. 
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concerns the isoscalar monopole components (C) and (D) in Fig. ([3]). The K"^ = 0^ 
of the giant quadrupole resonance (C) possesses also an appreciable monopole strength, 
therefore also the monopole resonance is split. The calculations for K"^ = 1^ are per- 
formed in such a way that the lowest solution is at zero energy. In that case all the 
spurious strength is concentrated there and the remaining states correspond to internal 
excitations of the nucleus. In Fig. Q the result of the electric dipole in ^'''^Yb is shown. 
Here a Gaussian with of 1.5 MeV FWHM was folded in all the levels in order to simulate 
the single particle widths. All results are in fair agreement with the experiments. 




Figure 4: Distribution of the = and 1 part of the B(E1) transition strength. Here 
a Gaussian with 1.5 MeV FWHM was folded into the discrete levels. 



5.4 Microscopic structure of the giant resonances 

In order to get more insight into the nature of the microscopically calculated collective 
states one may compare the microscopic transition densities with the classical picture. 
This quantity corresponds most closely to the density change of a classical vibration 
at maximum elongation. The transition density is defined in the intrinsic coordinate 
system as: 

P%,s{r = J2^Ur)xl'^M^) (39) 

\,fj, 
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where 92^(r) and x ^ire the single particle wave functions and QRPA amplitudes, 
respectively; K denotes the corresponding quantum number and S distinguishes the 
different solutions belonging to the same K. These transition densities have the same 
symmetries as Y\x and therefore only one quadrant of the z,x plane is drawn in the 
following examples. As a general result one finds that the transition densities of the 
giant resonances in different nuclei are very similar all over the deformed rare earth 
region. Therefore the result for the split giant quadrupole resonance in ^^"Yb which 
is shown in Fig. d5| is typical for all well deformed nuclei. The K'" = 2"'" excitation 





Figure 5: Microscopic transition densities between the ground state and the spit K- 
components of the giant quadrupole resonance in units of 10~^[/m"'^]. The 
dotted line denotes the boundary of the nucleus where the density attains half 
of his inside value. 
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corresponds closely to a classical 7- vibration and also for the K'^ = 1"^ one gets the 
expected pattern. The K'^ = 0"*" looks very similar to the classical /3-vibration. However, 
there is also change of the density inside the nucleus, which corresponds to a classical 
compression mode. As this compression vibration is predominately along the z-axis the 
authors called it axial breathing mode. The state "F" in Fig. [Hjwith K'^ =0"'' at 18 
MeV corresponds to the breathing mode in a spherical nucleus. As the corresponding 
compression vibration [H] is perpendicular to the symmetry axis, the authors called it 
radial breathing. This explains why one obtains a splitting of the breathing mode in 
deformed nuclei. 



6 Isomer Shifts 

An especially nice application of the second order response theory is the calculation 
of the change of nuclear charge radii due to rotation. These very small effects have 
been measured by two different methods: (I) applying the Mossbauer effect 01] and (II) 
using muonic atoms [45] . In both cases one observes the nuclear 2"^ — >■ O"*" rotational 
7-transitions in deformed even nuclei. From the view of the liquid-drop model the 
change of the radii has to be always positive due to the stretching effect. However the 
experiments showed positive and negative 5 (r^) in different nuclei, which ruled out this 
explanation. This paradox, which had been controversially discussed in the seventies, 
was solved with help of the extended Migdal theory described in Section 2.5. From the 
Mossbauer experiment the product 

^^Moes.6a.erOc5|*(0)|2<^(r2) (40) 

can be extracted, with 5|^'(0)|^ being the difference of the electron densities at the 
emitting and the absorbing nucleus and 5 (r^) the change of the mean-square charge 
radius 



^{^^) = \j d\rl5p{v) (41) 

where (5/?(r) denotes the change of the nuclear charge density upon excitation and Z is 
the charge. In the case of muonic atoms the shift is given as: 

5E^ = J d\ V^{r) 6p{r) (42) 

where is the Coulomb potential of the muon in the Is state. As one calculates the 
change of the density in the second order response one is able to calculate both shifts 
simultaneously. In muonic atoms one additional complication arises, because the 2^ 
state is split with a strong Ml transition between the two magnetic hyperfine doublets. 
However, the hyperfine splitting can be calculated within the same theory |36j. The 
authors [2T] calculated the excitation energy of the 2"'"-states, the change of the charge 
and mass radii and muonic isomerhifts. The results are in a fair agreement with the data. 
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Figure 6: The calculated change Sp^ of the occupation provabilities of the Nilsson levels 
near the Fermi energy for protons in "^^^Sm. The levels are drawn in a schematic 
way, equally spaced in the order of increasing energies e^. Open circles refer to 
N=5, full circles to N=4. The thin line gives the occupation probabilities in the 
ground state (scale is different from 6p^ 
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The most remarkable result, however, is their explanation of the physical mechanism 
which gives rise to positive and negative (5(r^). The effect is not connected with the 
collective motion but with the antipairing effect and the single particle structure near 
the Fermi edge. The antipairing effect tends to depopulate levels just above the Fermi 
edge in favor of levels below. It is important to mention that the single particle levels are 
spit due to the deformation. This has the consequence that near the Fermi edge one has 
proton states with the main quantum number = 4 and = 5 which have different 
radii. The change of the radii depends therefore only on a few levels at the Fermi edge. 
In Fig. (j?! the changes of the occupation probabilities of the Nilsson levels for ^^^Sm is 
plotted, where the change of the radius is positive. In ^^*^Dy the change of the radius 
is negative and the corresponding changes of the occupation probabilities are shown in 
Fig.([8|). In the case of "'^^■^Sm shown in Fig.(|8|), for three N=5 levels (with larger radii) 




Figure 7: Same as in Fig.(j6]), but for ^^°Dy. 

around the Fermi edge the occupation probabilities are increased and for two N=4 levels 
they are decreased which result in a positive 6 (r^). Depopulation of a larger N=5 level 
in favor of a smaller N=4 level leads to a negative isomer shift. This obviously occurs 
for ^^'^Dy due to the |~[523] proton level just above the Fermi energy. After this very 
simple explanation, the experimentalists lost the interest in these investigations. 
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7 Magnetic Excitations in Deformed Nuclei 



Magnetic excitations are calculated in the same way as the electric states discussed in 
Section 4. There is a large body of experimental data on Ml-transitions, that have been 
reviewed e.g. in [33 and an equally large number of theoretical investigations reviewed 
by Zawischa [IH]. There are two different modes: (I) one which is dominated by the 
orbital transitions and (II) spin-flip transitions which are known from spherical nuclei. 
While for the first class of transitions the experimental data and the theoretical results 
are well established, the interpretation of these states in terms of a collective model 
(scissor modes) has caused some discussions |48j. Like in the case of the electric /?- 
and 7-vibration of section 4, the low-lying states seem to have less resemblance with 
the collective model, compared with a predicted high lying, very collective resonance 
at around 23 MeV ^49j. For the spin-flip states such a problem did never appear. The 
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Figure 8: Results of a QRPA calculation of the B(Ml)t strength distribution in ^^^Gd 
with the full effective interaction and the bare magnetic operator. Spin- flip and 
orbital transitions are shown separately. 



18 



predominant orbital states are energetically lower compared with the spin-flip states. 
As a typical example, in Fig. ^ the results of a QRPA calculation for ^^^Gd are 
presented. Due to the mixture of spin-flip and orbital angular momentum states the 
Migdal parameter gQ and /q enter into the calculations. Both ph-force-components are 
repulsive so that the strength in both cases is shifted to higher energies. 




Figure 9: Spin-flip strength distribution in ^^"'Sni. The experimental data[52] are com- 
pared with QRPA results. Gaussian have been folded in the theoretical results to 
produce a continuous distribution. The total theoretical result is given by the full 
curve. The broken curve in the upper part represents the incoherent sum of the 
proton and neutron contributions shown in the lower part. 

There is a clear separation between the orbital components, the so called scissors- 
modes and the energetically higher spin- flip transitions. In Fig. ^ a comparison of 
theoretical results with the data for ""^^^Sm are given [5Uj. This double bump structure 
is characteristic for all well deformed rare earth and actinide nuclei |,51j . 
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8 Derivation of the Landau parameters 



The Landau-Migdal parameters can be derived from an underlying many-body theory 
of nuclear matter and finite nuclei. One has to start from a theory of the nuclear ground 
state and perform a functional differentiation of the ground state energy of nuclear 
matter with respect to the quasi particle occupation numbers, as was suggested first by 
G. E. Brown and S. O. Backmann |53 t I54 t 155] . Several deep insights into nuclear physics 
have been obtained this way. First of all, it turned out that the spin-isospin dependent 
part of the interaction, the famous parameter could be understood quantitatively from 
the one-pion exchange and from a reasonable assumption about short-range correlations. 
This is ultimately due to the long-range character of the pion-exchange which allows to 
treat the effects of short-range correlations in a simplified way by essentially removing 
the short range part of the pion-exchange. 

An important phenomenological generalization of the Landau-Migdal interaction has 
resulted from this observation j56| 157] . The Stony Brook Juelich ansatz augments the 
conventional Landau-Migdal parametrization by the explicit one-pion and one-rho ex- 
change and Vp which are folded by a correlation function J7c(q): 

G'{q) =J 7^ [^-(k) + Vp{k)] 17e(q - k) + 6G',a ■ a'r ■ r' (43) 



with 
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Oe(q) = (2vr)^5(q) - — 5(|q| - g,). (44) 

Here qc = 3.93 denotes the inverse of the Compton wavelength of the uj meson. A 
parameter 6G'q is introduced to account for a small correction to the Landau-Migal 
parameter G'{q) that is not produced by the explicit meson-dynamics. 



In the upper part of Fig. (10) the q-dependence of G' is plotted. One realizes that 
at small momentum transfers the central part of the spin-isospin interaction is strongly 
repulsive and for larger momentum transfers it is weak. The central part of the vr - 
meson and p -meson contributions have the same sign, whereas the tensor parts have 
the opposite signs. The p -meson exchange therefore acts as a natural cut-off for the 
strong tensor component of the one-pion exchange. All details of the calculations are 
given in the original publications [56] . 

The unnatural parity 12~ and 14~ magnetic high spin states discovered experimentally 
by Jochen Heisenberg and his collaborators at the BATES electron scattering facility 
lool m are a striking example to illustrate the momentum dependence of the 

generalized Landau-Migdal interaction. The cross sections peak around g ~ 2 [fm~^], 



see Fig. (11). The energies of the two 12~ and the one 14~ states are close to the 
experimental ph-energies but the cross sections are only half of the shell model prediction. 
The explanation for this surprising result was given in Ref. [SUj. As the spin-isospin 
interaction in the relevant momentum range is essentially zero, the RPA solutions are 
close to the unperturbed ph-energies. In their extended model the authors included 



20 



0,8 
0,6 
0,4 
0,2 
0,0 









Central part 





















0,0 
-0,2 
-0,4 
-0,6 











.... 




Tensor part 







q[fm-^] 3 



Figure 10: Momentum dependence of the spin-isospin parameter G{q{ in units of 
[SOOMeT^/m^]. The full line is the complete model of Eq. the dashed 

hne in the lower part is the correlated one-pion exchange only. 
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Figure 11: Inelastic electron scattering cross sections at O = 90deg and 9 ~ 160 deg of 
the magnetic high spin states in ^°^Pb. Experiments are compared with the 
RPA results (dashed lines) and the extended model. Calculations have been 
done in DWBA. 
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also the effects of the low-lying phonons within the so called core coupling RPA which 
provides an explanation for the reduction of the cross sections. 

While magnetic modes with high angular momenta do not show collectivity, magnetic 
modes of small multipolarities may show surprisingly large cross sections. The most 
outstanding mode of this kind is the Gamov- Teller resonance discovered by Charles 
Goodman and Dan Horen at the Indiana cyclotron facility in charge-exchange reactions 
[SniEI]. In Tab. ([T]), the averaged ph-energies eph are compared with the RPA excitation 
energies E^p^ for the 0^, 1^ 1^ and 2~ spin-isospin modes in ^"^Pb. In the first three 
cases, the energy shift AE is of the order of 5 MeV. The 2~ result is qualitatively 
different, however. Here the ph- force is weak as the transition density is peaked at 
larger momentum transfer, and as a consequence one obtains four states which are only 
little shifted from the uncorrelated ph-energies. The comparison with the data shows a 
fair agreement as far as the mean energies are concerned. 





eph [MeV] 


Erpa [MeV] 


l^E [MeV] 


Ee^p [MeV] 




21.8 


26.6 


4.8 


25.1 ± 1.0 


1+ 


13.1 


18.9 


5.8 


19.2 ± 0.2 




21.5 


26.3 


4.8 


25.1 ± 1.0 




20.7 


22.6 


1.9 




2- 


20.1 


21.2 


1.1 


25.1 ± 1.0 




23.3 


24.5 


1.2 






28.0 


28.8 


0.8 





Table 1: Charge-exchange resonances in ^'^^Pb. 



The charge exchange resonances have relatively large widths which can not be obtained 
within a Iplh approach but one has to include higher configurations. In Fig. (12) 
an example for such more involved calculations is shown. The authors f57| extended 
the conventional RPA approach and include Iplh as well as 2p2h-configurations in a 
consistent way. This calculation not only reproduced the known experimental energy 
and width but it also predicted a long tail up the 50 MeV, where nearly half of the 
strength is hidden. A detailed discussion of spin-isosopin modes is given in Ref. [03] . 

In the early 1970's, there were speculations about the existence of a pion condensate 
in nuclei. At large nuclear densities, the spins of the pions were supposed to align as 
a consequence of the one-pion exchange interaction. The investigation of high energy 
heavy ion collisions were suggested as an experimental tool to study the supposed phase 
transitions, despite warnings that the time-scale for performing a phase transition does 
not match with the time two high energy ions need to pass each other [M]. More over, 
one had to realize that the pion-exchange is always accompanied by short-range nuclear 
correlations which generate a Landau- Migdal parameter g^, and thus suppress the onset 
of a phase transition [64^ . 

In the case of the Landau parameter /q, microscopic derivations had to face severe 
challenges. Landau had related the parameter /o to the compression modulus of nuclear 
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Figure 12: Gamow- Teller resonance in ^°*Pb calculated in an extended RPA model [57] , 
The arrow indicated the experimental mean energy. 



matter, 

i^ = 3^(l + /o). 
m 

Inserting the empirical value K = 210 it 30MeV, one finds that /o cannot be negative. 
Moreover, the isotope shifts in heavy nuclei clearly rule out negative values of /o P] • 

Brueckner's theory of nuclear matter produces values of the spin- and isospin- inde- 
pendent parameter violating Landau's stability criterion /o > —1, see Table 

This is an expression of the fact that Brueckner theory does not saturate nuclear 
matter at the empirical saturation point, but at much larger densities. The advent of 
Walecka's relativistic mean field theory suggested a new mechanism to saturate nuclear 
matter at the correct saturation point [U^ . Indeed, incorporating the relativistic effects 



Table 2: Landau parameters of nuclear matter for a Fermi momentum of kp = 1.36fm~^, 
obtained by using the one-boson exchange potential HEA 152J as the two-nucleon 
interaction and different versions of Brueckner theory to derive the binding energy 
of nuclear matter. First row: standard Brueckner Hartee-Fock. Second row: 
The Dirac spinors of the bare two-nucleon interaction are modified Third row: 
the empirical parameters determined from linear response of finite nuclei. All 
parameters are given in units of 300 McV fm"^. 



method 


/o 


/o 


90 




BHF, non- relativistic 


-1.17 


0.32 


0.2 


0.63 


BHF, relativistic 


-0.72 


0.42 


0.12 


0.63 


empirical 


-h0±0.2 


0.8 


0.2 


0.9 
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suggested by Walecka into a Brueckner calculation, one is able to produce a Landau- 
Migdal parameter /o which signals stability of the nuclear ground state, but quantita- 
tively, it does not agree with the empirical one. This clearly shows that the Walecka 
approach is incomplete [15]. A hint to the relevance of three-body interactions came 
from phenomenology. The effective interactions of the Skyrme type rely on medium- 
dependent effective interactions which may be interpreted as effects due to three-body 
forces. The different phenomenological Skyrme-like parameterizations allow a dramatic 
variation in the resulting Landau-Migdal parameters, however, which blocked progress. 
Indeed, instead of deriving the Landau parameters from an underlying many body ap- 
proach, one rather used the empirical Landau-Migdal parameters to constrain the free- 
dom available in the choice of appropriate parameterizations of the effective medium 
dependent interactions [BS]. For progress, a theory of the three-nucleon interactions was 
required. In the 1980's the then new field of hadron physics emerged and improved our 
understanding of hadronic interactions. Effective Field Theories for pion-pion scatter- 
ing, pion-nucleon interactions, and eventually nucleon-nucleon scattering were developed 
and provided a systematic expansion scheme which restricts the vast number of three- 
nucleon interactions allowed in older theories. A recent summary is given in Ref. [67]. 
In next-to-next-to-leading order, there are exactly three types of three-nucleon interac- 
tions. First theoretical investigations of nuclear matter based on those interactions have 
become available and are expected to make a major impact on the theoretical command 
of the Landau-Migdal parameters [681 ESj • A successful description of neutron matter 
based on chiral three-nucleon interactions has been published and promises a solid basis 
for the theory of neutron-rich nuclei |7U] . 



9 Conclusion 

We have reviewed the application of the conventional TFFT to deformed nuclei in the 
rare earth and actinide region. The electric as well as magnetic states have been cal- 
culated. The quasi particle in the definition of Landau are the Nilsson single particle 
levels. Whereas in spherical nuclei the corresponding single particle levels are 2j -|- 1 
times degenerate, this degeneracy is lifted due to the deformation. For that reason the 
density of single particle states is much higher and one does not need e.g. an additional 
splitting due to low-lying collective states. The theoretical results are in general in 
good agreement with the experimental data. The low-lying (3- and 7- vibrations are not 
very collective, therefore the theoretical results depend strongly on the single particle 
energies. We also reviewed an extension of Migdal's theory, the so called second order 
response theory. An especially nice example are the isomer shifts of rotational states 
which we discussed in some details. Finally we discussed more recent developments in 
the frame work of effective theories. 



25 



10 Acknowledgement 

It is a pleasure to thank S. Drozdz, V. Klemt, J. Meyer-ter-Vehn, P. Ring, G. Tertychny"'", 
J. Wanibach and D. Zawischa for the collaboration in developing and extending the 
Landau-Migdal theory in the past, S. Avdeenkov, S. Kamerdziev, N. Lyutorovich and 
V. Tselyaev for their collaboration up to present. We thank the Deutsche Forschungsge- 
nieinschaft for the grants DFG: 436 RUS 113/994/0-1 and DFG: 436 RUS 113/806/0-1. 
JS thanks the Foundation for Polish Science for financial support through the Alexander 
von Humboldt Honorary Research Fellowship. 

References 

[1] A.B. Migdal, Theory of Finite Fermi Systems and Application to Atomic Nuclei, 
Wiley, New York, 1967. 

[2] J. Speth, E. Werner and W. Wild, Physics Reports 33, (1977) 127. 

[3] J. Speth Z. Physik 239 (1970) 249. 

[4] P. Ring and J. Speth, Nucl. Phys. A235, (1974) 315. 

[5] V. A. Khodel and E. E. Saperstein Physics Reports 92, (1982) 183. 

[6] S. Droidi, S. Nishizaki, J. Speth and J. Wambach, Physics Reports 197, (1990) 1. 

[7] S. Kamerdzhiev, J. Speth and G. Tertychny, Phys. Rep. 393, (2004) 1. 

[8] V.I. Tselyaev, Sov.J.Nucl.Phys. 50, (1989) 780. 

[9] V.I. Tselyaev, Phys. Rev. C75, (2007) 024306. 

[10] F. Griimmer and J. Speth J. Phys. G: Nucl. Part. Phys. 32 (2006) R192. 

[11] S. Krewald and J. Speth, Int. Journal of Mod. Phys.18 (2009) 1425. 

[12] V.I. Tselyaev, J. Speth, F. Griimmer, S. Krewald, A. Avdeenkov, E.V. Litvinova 
and G. Tertychny, Phys. Lett. B653, (2007) 196. 

[13] D.Vretenar, N. Paar, P. Ring and G.A. Lalazissis, Nucl.Phys. A692, (2004) 281. 

[14] E.V. Litvinova and V.I. Tselyaev, Phys. Rev. C75, (2007) 054318. 

[15] D.Vretenar, T. Niksic, N. Paar and P. Ring, Nucl.Phys. A731, (2001) 496. 

[16] N.Ryezayeva et al. Phys. Rev. Lett. 89, (2002) 2727502. 

[17] D. Vretenar et al.Phys. Rev. C65, (2002) 021301. 

[18] S. Krewald, K. Nakayama and J. Speth, Phys.Rep. 161 (1988) 103. 



26 



[19] L. P. Kadanoff and G. Baym, Quantum statistical mechanics, W. A. Benjamin, Inc, 
New York, 1962. 

[20] W. Brenig and H. Wagner, Z. Physik 173 (1963) 484. 

[21] J. Meyer and J. Speth, Nucl. Phys. A203 (1973) 17. 

[22] M. G. Urin and D. Z. Zaretsky, Nucl.Phys. 35 (1961) 219. 

[23] M. G. Urin and D. Z. Zaretsky, Nucl.Phys. 47 (1963) 97. 

[24] B. L. Birbrair, Nucl.Phys. A212 27. 

[25] L. D. Landau, JETP 3 (1957) 920; 5 (1957) 101; 8 (1959) 70. 

[26] A. L. Larkin and A. B. Migdal, JETP(Sov.Phys.)17 (1963) 1146. 

[27] B. L. Birbrair, Nucl. Phys. A108 (1968) 449. 

[28] S. P. Kamerdzhiev, Sov. J. Nucl. Phys. 9 (1969) 190. 

[29] J. Speth, Nucl. Phys. A135 (1969) 445. 

[30] L. P. Gorkov JETP(Sov.Phys.)7 (1958) 505. 

[31] N. N. Bogolyubov, ZhETF(UDSSR) 67 (1958) 58. 

[32] M. Baranger, Phys.Rev.120 (1960) 957. 

[33] S. T. Belyaev, Nucl. Phys. 64 (1965) 17. 

[34] J.Meyer and J. Speth, Phys. Lett. B39 (1972) 330. 

[35] A. B. Migdal, Nucl. Phys. 13 (1959) 655. 

[36] B. L. Birbrair, Nucl. Phys. A108 (1968) 449. 

[37] J. Meyer, J. Speth and J. H. Vogeler, Nucl. Phys. A193 1972 60. 

[38] A. Bohr and B. Mottelson, Nuclear Structure, Vol.11 (Benjamin, reading, Mass., 
1975) 

[39] J. H. Vogeler, Nucl.Phys. A133 (1969) 289. 

[40] C. Gustafson, I. L. Lamm, B. Nilsson and S. G. Nilsson, Ark. Fys. 36 (1971) 69. 
[41] D. Zawischa, J. Speth and D. Pal, Nucl.Phys. A311 (1978) 445. 
[42] E. G. Fuller and E. Hayward, Nucl.Phys. 30 (1962) 613. 

[43] M. N. Harakeh and A. van der Woude, Giant Resonances, Oxford University 
Press, Oxford, 2001. 



27 



[44] G. K. Shenoy and F. W. Wagner, Mossbauer Isomer Shifts North-Holland, Ams- 
terdam, 1978. 

[45] G. Goldring and R. Kalish, Hyperfine Interaction in Excited Nuclei, Gordon and 
Breach, New York, 1971. 

[46] J. Meyer, P. Ring and J. Speth, Phys. Rev. C7 (1973) 1803. 

[47] K. Heyse, P. v. Neumann-cosel and A. Richter, Rev. Mod. Phys. 82 (2010) 2365 

[48] D. Zawischa, J.Phys. G:Nucl. Part. Phys. 24 (1998) 683. 

[49] D. Zawischa and J. Speth, Phys. Lett. 219B (1989) 529. 

[50] D. Zawischa and J. Speth, Nucl.Phys. A569 (1994) 343c. 

[51] D. Zawischa, M. Macfarlane and J. Speth, Phys. Rev. C42 (1990) 1461. 

[52] D. Prekers et al., Phys.Lett. B244 (1990) 178. 

[53] G.E. Brown, Rev. Mod. Phys. 43 (1971) 1. 

[54] S.O. Backmann, Nucl. Phys. A120 (1968) 593. 

[55] S.O. Backman, A. D. Jackson and J. Speth, Phys.Lett. B56 (1975) 209. 

[56] J. Speth, W. Klemt, J. Wambach and G. E. Brown, Nucl.Phys. A343 (1980) 382. 

[57] S. Drozdz, V. Klcmt, J. Speth and J. Wambach Phys.Lett. 166B (1986) 18. 

[58] J. Lichtenstadt, J. Heisenberg, C.N.Papanicolas, C.P. Sarget, A.N. Courtemanche 
and J.S. McCarty, Phys.Rev.Lett.40 (1978) 1127. 

[59] S. Krewald and J. Speth, Phys.Rev.Lett. 45 (1980) 458. 

[60] C. D. Goodman et al. Phys.Rev.Lett.44 (1980) 1755. 

[61] W.G.Love, in The (p,n)Reaction and the Nucleon-Nucleon Force, edited by C.D. 
Goodman et al. Plenum, New York, 1980, p. 30; F. Petrovich, ibid., p. 135. 

[62] K. Holinde, K. Erkelenz, and R. Alzetta, Nucl. Phys. A198 (1972) 598. 

[63] F. Osterfeld, Rev. Mod. Phys. 64 (1992) 491. 

[64] S. Krewald and J. W. Negele, Phys.Rev.C21 (1980) 2385. 

[65] B.D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16 (1986) 1. 

[66] K. F. Liu, H. Luo, Z. Ma, Q. Shen, Nucl. Phys. A534 (1991) 1. 

[67] E. Epelbaum, H.-W. Hammer, and Ulf-G. Meifiner, Rev. Mod. Phys. 81(2009) 1773. 



28 



[68] S. K. Bogner, R. J. Furnstahl, A. Schwenk, Prog. Part. Nucl. Phys. 65 (2010) 94. 

[69] J. W. Holt, G.E. Brown, J. D. Holt, and T. T. S. Kuo, Nucl. Phys. A785 (2007) 
322. 

[70] K. Hebeler and A. Schwenk, Phys. Rev. C82 (2010) 014314. 



29 



